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abstract 


The J-integral for q,ua si -iso tropic ( 0 /+ 4^/90)23 
crossplied (0/90), laminates has been evaluated using 
finite element technique. Constant strain traingular finite 
elements have been used. The program gives displacement 
at all the nodes and stresses and strains in all the elements. 
Stresses are obtained for all the plies and J-integral can 
be obtained for all the laminae as well as for the laminaste. 
Three crack configurations namely, single edge 
notched (SW), double edge notched (DSF)and centrally 
located crack specimens have been analysed,' 

For the formulation of the problem, an equivalent 
stress-strain matrix for the laminate has been defined such 
as to model the laminate by a homogenous and orthotropic 
material, A stepwise method has been used to predict 
laminate failure. To account for effect of stacking 
sequence, models with no load relaxation and complete load 
relaxation in failed plies have been investigated. 

The results obtained by present study compare quite well 
>/ith the experimental investigation on SEN specimens. 
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INTRODUCTION 

1.1 IIWRODT'CTION TO FRilCTURE OF COMPOSITE MATERIALS 

Composite materials are engineering combination of 
two or more materials to achieve certain physical properties 
not realizable by constituent materials individually. 

Composite materials have high strength and stiffness 
coupled with low density. In addition, they have good 
vibration, fatigue and corrosion resistance, low heat 
conducjiivity, good electrj cal insulation and favourable 
cost effectiveness. Manufacture of composites requires 
less labour and generates less waste. It has the 
additional advantage of easy processibility. 

Since the composite materials are being increasingly 
used in many sophisticated applications, material 
selection criteria and design procedures are required. 

Global properties like deformation and modulus are well 
understood. Strength criteria, although not well ; 

established, are generally agreed upon. So, research 
emphasis is towards understanding of fatigue and fracture 
response of composites. 

Fracture mechanics has now been widely accepted 

as a useful discipline, which deals with inherent macroscopic i 

[ 

cracklike defects in materials. It deals with the study : 

of fracture of materials by crack initiation and propagation J 
under various service conditions. It aims at providing i 

I 

guidelines for tolerable size of cracklike defects and | 

\ 

[ 
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inspection requirements for avoiding catastrophic failures. 
Fo 3-> metals and ceramics, developments in fracture mechanics 
has led to establishing methods of predicting fracture 
toughness. The fibre reinforced composites also seek 
reliable means of characterizing fracture behaviour as a 
pa.rt of establishing the fundamental properties. 

Due to heterogeneity of the composite materials 
and complexity of interaction between constituents , 
fracture is accompanied by several processes. The. sc 
processes include fibre brealiage, fibre extraction, 
debonding between fibres and matrix, delamination, fracture 
of matrix etc. In some- cases composites fracture during 
stretching along the fibres by breakdown of integrity. 

Even when steady growth of crack is observed, they differ 
substantially from the cracks in traditional materials. 

The problem of characterizing fracture properties of 
composite materials is challenging as the crack propagation 
in them may not be self -similar, D”c to these complexities 
parallel or combined experimental and theoretical analyses 
is required. 

1*2 LITERATURE REVIBV/ 

In fracture mechanics two basic approaches, namely, 
the energy and stress approaches have emerged. Thu energy 
approach is due to Griffith [l], Hg derived a fraeture 
criterion by taking into account total change in energy 
of a cracked body. The stress approach is due to Irwin [2], 
who relied on stress analysis of a cracked body. 
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Definitive features of fracture in homogeneous 
isotropic solids vrere identified by Griffith [l] and 
extended by Orowan [3] and Irwin [2] for engineering 
materials, Irwin carried out stress analysis by assuming 
the material to be linear elastic. This approach is known 
as linear elastic fracture mechanics (LEFM) . 

Rice [4,5] used a different energy approach 
applicable to linear as well as nonlinear elastic materials 
and under conditions of large scale yielding also. He 
used a two dimensional energy line integral around crack 
tip, now better knoim as J-integral, This integral 
characterizes the crack tip. Basis of J-integral is 
provided by the works of Hutchinson [6] and Rice and 
Rosengren [7] .. I'fcClintock [3] has demonstrated that 
crack tip stress and strain tlelds can be described in 
terms of J-integral, The use of J-integral as an elastic- 
plas-uic fracture criterion has been discussed by Brobcrg[9]* 
Bagley and Landes [IO-I2] through extensive experiments 
have shown applicability of J-intcgral as fracture 
criterion for metals. 

Several researchers applied linear elastic frSLcture 
mechanics to composite material and found that fracture 
toughness was not independent of initial crack length and 
they used a two parameter LEFM model. Hr a consideration 
of hole size effect. Fuismer and Whitney [13,14] proposed 
alternative approaches for fracture of laminated composites# 
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TI:ese two approaches are point stress criterion and average 
stress criterion, Giare [15] and Agarwal and Giare [16] 
found notched strength of short fibre composites to be in 
agreement with Whitney-Muismer criteria. 

The J -integral approach was attempted on glass fabric 
reinforced polyester resins by Smith, Green and Bowyer [l/] . 
It has been used to short glass fibre composites by Agar'.Aral, 
Patro and i^mar [18,19], Srinivasan [20] has used 
J-integral for quasi -isotropic and crossplied laminates. 

Babu [21] in a theoretical study has evaluated 
J-integral for short fibre and unidirectional composites 
using finite element method. He used a 2-dimensional 
displacement finite elomont analysis using constant strain 
traingular elements. J-integral for randomly oriented 
short fibre composites and unidirectional composites was 
evaluated for piano strain and plane stress cases. 

He conducted the analysis for three different crack 
configuration, namely, single edge notched, double edge 
notched and contrally located crack. His check for path 
independency of J-integral was satisfactory, J for short 
fibre composites are in good agreement with experimental 
works of Patro [22]. 

1,3 gGOPE OF PEESENT VDRIC 

The present studies arc an attempt to predict the 
fracture behaviour of composite laminates using J-integral 
as a fracture criterion. Symmetric laminates have boon 
analysed using equivalent stiffness matrix for the laminates. 
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J -integral has been calculated for individual laminae 
through lamina stiffness matrix. Attempt has been made to 
predict first ply failure and analysis after first ply 
failure has been carried out by modified stiffness matrix. 

The analysis has been carried out for quasi -isotropic 
[0/± 45 / 90 ] 2 g crossplied [ 0 / 90 ]]^g laminates for three 
different crack configurations, namely, single edge 
notched, double edge notched and centre located cracked.' 

Formulation of the problem is given in chapter 2. 
Finite element solution technique adopted for the present 
study is also discussed in chapter 2 . Chapter 3 deals 
with the numerical results obtained. The various models 
adopted to predict fracture behaviour of laminates are 
also discussed in this chapter, in the light of experimental 
fifidings. Conclusions ana some suggestion for future work 
are included in chapter 4* 3Iho computer program used for 
the present study has been given as an appendix.' 



CHAPTER II 


J -INTEGRAL AS FRACTURE CRITERION FOR COMPOSITE mTERlALS 


2. 1 II'ZCRODUCTION 


Most of the research work on fracture of composite 
materials j so far uses linear elastic fracture mechanics 
(LEFM) concepts and requires accurate stress analysis at 
the crack tip. Due to stress singularity at crack tip, 
this is difficult even for homogenous isotropic materials* 
In heterogenous composites, the situation becomes more 
difficult due to complex crack tip geometry. 

The J-integral is a parameter which can be evaluated, 
due to its path independency, without focussing attention 
at the crack tip, but can still characterize the crack tip. 
J-integral does not require a prescribed crack extension 
trajectory and is hence more adaptable to characterize 
cracks in composite material. 

The J-integral, as proposed by Rice [4,5] is a two 
dimensional energy line integral defined as 

(Udy “T^ - - ■’ ds ) (2.1 ) 

where p is any contour traversing in counter clockwise 

sense and enclosing tho crack tip as shown in Pig, 2.1^ 

U the strain energy density is defined as 

e 

/> mn 


U = U ( ^ ) 

^ mn' 


^ ®iD 


Where c . . and are stress and strain tensors 

13 13 
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respectively* T^. is the traction vector and is given hy 
T . = e . . n . 

¥hcrc n^ is the normal vector to p . is the displacement 
vector and s the arc length along p . 

2.2 FOHiULATION OF PROBLEli 

For evaluation of -J -integral, stresses, strains and 
derivatives of displacement arc to b. found out* To achieve 
this end total potential energy stored in the laminate has 
been minimized by Raleigh -Ritz method [23]. The finite 
element program, then solves the problem mth prescribed 
boundary conditions to yield displiccracnt at all the nodes, 
their derivatives, strains and stresses* 

2.2.1 Expression for Potential Energy 

For symmetric laminates subjected to inplanc loads 
only, displacement and strain fields are same throughout 
the tld-ckness of the laminate, and thus donot vary from 
lamdna to lamina [24-]. 

Due to different fibre dirt ction , stress-strain 
matrix vary from lamina to lamina and honce the stresses 
in them arc not same oven with same strain field. 

Now if the stress-strain relationship of i^ 
lamina in a laminate be [D]j^ and its thickness bo t^, 
strain energy stored in it is given by 

i i ^ 
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where LeJ is the row matrix ®2 '^12J 2 i,nd 

is the column matrix. So the total strain energy stored in 

the laminate will he 

U = Z u =1 t.ej (^[ ^li (2.3) 

i i 


Now if the laminate is a flat plate 

dVj_ = dA t^ (2.4') 

and eqn, 2.3 can he written as 


U = 




(2.?) 


Now an equivalent stress strain matrix for the 
laminate is defined such that 


[D]e t 


c 


£ [D]i q 

i ^ 


where [d] is the equivalent stress-strain matrix and 
o 

is the total laminate thiclmess. 

Equation 2.5 can he rewritten as 


( 2 . 6 ) 


U = l_ej [d]^ dA 


(2.7) 


Now in the absence of body forces expression for the potential 
energy of the laminate can he written as 


I = 


L®J P ], 



t 

c 



I u as 
a 3 


( 2 . 8 ) 
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2 • 2 ♦ 2 Boun dazy C on di ti on 

Three different crack configurations namely single 
edge notch*' d, double odge notched, and centrally located 
crack plate.' have been analysed. Those plates have hoen 
subject.',! to mode I loading condition as shown in Fig. 2.2. 
Boundary condition in all of them are 


Alone AB 


y = 0, r =0 
> sy 


Along BC and DA 0=0; r = 0 (Traction free) 

Along CD 7 = 6 (applied displacement) 

and 'tr =0 
xy 

crack surfaces are stress free, 

2.3 FINITE ELEMENT SOLUTION TEGHMIQUE 

Expression for the potential energy of the laminate, 
Eq, 2.8, contains only first derivative of displacements. 
SO compatibility of u,v and completeness of 


u, V, 


Du 

ax 


B V 
ax ’ 


0 U 

dJ ’ 


are required [ 23 I. 


ay 


Traingular elements with throe nodes with following 
interpolation function satisfy the- above requirements. 

e f 

{u] = a +■ bz + oy = UJ 3 _ Mg WjJ I Uj 


Tl-‘us 


u- 


* V 0 

|vj = a* + b*x + c'y = [jlTj 


Np N^ J ■\ Vo \ 


j-1) 

2 
V. 


(2.9) 


where [^Nj^ N^ N^J are interpolation functions. 

Substituting equation 2.9 iu 2.8 and applying the Raleigh- 
Ritz method [ 23 ], one gets 



OCV' 



THREE CRACK 
WITH BOUNDAR 

I 





1 2 


[k] 


d' 


>(ne) 


= |p] 


(ne) 


( 2 . 10 ) 


e f f T 

where [k] = JJ [b] [d] [b] t clA = stiffness matrix. 


A(e) 


[B], 


c)Nj 
2 X 

0 


1 y 


0 




c- y 


;>iL j 

X i 


,(ne) 


and > d\ 


I / 

iJ 


i is the node number. 

So for traingular elements 


[B] = 


ah. 


ax 

0 


BN. 


1. 


^ y 


ne 


and i d 


3 nr. 


BN 


1 


ay 


an. 


Sx 


0 


aN. 


dx 


■^1 
y ^2 


u- 






0 


a N. 


a y 




B X 


0 




ay Bx 3 y 


0 


aF. 


a y 


a X 


( 2 . 11 ) 


( 2 . 12 ). 
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( 

I 


P 




P 




ds = Load matrix 


( 2 . 43 ) 


Ng = interpolation function along the houndaries of the 
finite element, 

Now the specimen to he analysed is divided into a 
finite number of traingular elements. 1 moderately fine 
element grid is maintained near crack tip. The size of 
elements is increased as it is moved away from crack tip 
because of less variation in stress. Element stiffness 
and load matrices are obtained and assembled. After 
applying appropriate boundary condition! displacement, 
derivatives of displacement and strains are obtained. After 
this using the lamina properties stresses in different 
laminae and consequently J is determined for all laminae 
separately. In view of the path independency, contour 
is chosen along the boundary of the specimen as it is 
convenient. 


2, if A CORRECTION FACTOH 

For homogenous and orthotropic material , relation 
between I and stress inteUvSity factor K is given by [25] 





(2,llf) 


t' 


where 



^11 ^22 ^ 1/2 j - ( 

2 


^22 

ail 


^ 1/2 ^ 


^66 ^ ^^12 


2a. 


11 


( 2 , 15 ) 



Ilf 


a. . are terms of the compliance matrix of the material. 

J- J 


how if the elastic modulus in loading direction he 


E 


E 


then J = ( — Y ) 


'1 


E 


E. 


how from the finite clement formulation only 

E 


factor is evaluated, so a correction factor 
introduced to get J value for tho laminate. 


E 


U . 

"Y IS 


( 2 . 16 ) 


a 



CHA.PTER III 


RESULTS AND DISCUSSION 

This chapter deals with the numerical results obtained 
for q.uasi-isotropic and crossplied laminates and their 
comparison with experimentally obtained values. An attempt 
has been made to predict laminate failure by using properties 
of constituent laminae. A stepwise method has been used to 
predict laminate failure. 

3,1 PROCEDURE FOR LAKENATE ANALYSIS 

For a given specimen, J-integral was evaluated at 
different displacements for the laminate as well as for the 
laminae j constituting the laminate. At small displacement 
none of the plies fails, A^ the load increases, plies start 
failing. A ,-ily was assumed to fail when the applied 
displacement exceeded a critical displacenent for the ply. 

The final laminate fracture is assumed at the critical 
displacement for the laminate obtained exp crimen tsOIy, Tie 
basic lamina properties and the critical displacements for the 
purpose of present investigation wore obtained experimentally 
and are given in Tabic 3»1, 

After a ply fails, its contribution to the overall 
laminate behaviour changes, Tliis change is incorporated in 
the analysis through modification of its stiffness matrix. 
There are different approaches which can be used to modify 
the stiffness matrix [2.4], One approach suggests that the 
stiffness matrix of a failed lamina should bo made zero. This 
would mean that the lamina would carry no additional load in 
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any direction. This will be a conservative approach since, 
for example, when a ply fails by transverse tensile failare, 
it may not carry any additional load in transverse direction 
but can certainly carry additional load in the fiber direction 
A more practical approach under such circumstances would be 
to assume that while thu transv-rsc tonsils modulus and shear 
modulus of a failed ply arc zero but the longitudinal modulus 
remains unaffected, Th*, latter approach has been used in the 
present analysis. 


TABLE 3^1 

LAI’CCNA ELASTIC CONSTANTS 


Longitudinal modulus, Ej^ 


31.9 

Transverse- modulus, 

— 

9.7 GPa 

Major poisson's ratio, V 


0.25 

Shear modulus, 


3.2 GPa 

Lamina.te s 

Critical displacement 
(mm) 

0° unidirectional 


0.88 

90 ° unidirectional 


0,08 

4-5° unidirectional 


0.10 

Quasi -iso tropic, (0/+ 45 / 90)23 


0 , 75 * 

Crossplied, (0/90)]^ ^ 


0.375* 


♦ From reference 20 
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The J-integral evaluated through a step-vTise procedure 
may he represented hy a piecewise continuous curve as shown 
in Fig, 3.1, v/here each discontinuity point represents a 
ply failure. In the segment heti'/een two discontinuities the 
J-integral, (J)^^, can be written as follows 

(J)^ = ^ l)j_ (3.1) 

Where (J)^_p is the value of J at the (i-1)'^^ 
discontinuity and j)^ is the increamental J in i^^ segment. 
When all the plies are intact, the J-integral is evaluated 
using the original stiffness matrix for the laminate. The 
increamental J is calculated using modified stiffness 
matrix. 

Another point which must be considered at this point 
is what happens to the load carried by the ply when it falls 
That is, when a ply fails what is the extent to which the 
load carried by it drops. In a real situation, it may be 
between the two extreme situations, namely, one in which the 
load in felled ply drops to zero ( complete load relaxation) 
and two, the load doesnot drop at all ( no load relaxation). 

It can be easily shown that the later condition would 
show greater increament of J compared to the former condition 
for the same displacement increament. In the present 
investigations, results have been obtained for both the 
conditions, 

3, 2 J-INTEGRAL FOR LAMINATES 

J-integral ms evaluated as a function of applied 
displacement (&) for q.uasi-isotropic and crossplied 
laminates. The ratio of crack length to specimen width was 
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FIG. 3.1 J vs 6 PLOT FOR - A LAMINATE 
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yarried from 0,1 to 0,7. The following crach configurations 
were used: 

i) Single edge notched (SEN) specimens 

ii) Double edge notched (DEN) specimens 

iii) Centrally located notch specimen 
J YS S curves for SEN specimens are shown in Figs. 3.2 and 
3.3 for quasi-iso tropic and crossplied laminates respectively, 

J is plotted against in Figs, 3.4- and 3«5. Since J vs 
plots are straight lines passing through the origin, it is clear 
that J-dntegral is proportional to the square of -applied 
displacement for all these cases. In fact, this behaviour is 
expected. In the present analysis, it has been assumed that 
the laminate stress-strain ( or load-displacement) behaviour is 
linear. The load -displacement cuinres for a laminate with 
different crack lengths may be represented as shown in 
Fig, 3*6. For the linear relationship, the potential energy 
stored is given by 

U = I c (3.2) 

where ’c’ is the slope of cruve. It should be 
noted thn.t c is a function of crack length and type of 
crack, but is independent of applied displacement ’ ® * , 
Therefore I-integral at applied displacement * & ’ will be 

j= - i = ^£2 ,{3^3) 

da ^ 3 a 

where k = - ^ is again a function of crack configuration 

zd a 

only. 



, KJ/m^ 



DISPLACEMEN 


FIG. 3.2 J vs 6 PLOT FOR 


KJ/m^ 



DISPLACEMENT , mm | 

FIG. 3.3 J vs 6 PLOT FOR DIFFERENT (a/w) | 



KJ/m 





, KJ/m 



(DISPLACEMENT)^ , mm^ 
FIG. 3.5 J vs 6^ PLOT 
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I - Crack Length 'a^ 
II - Crack Length 
02 



DISPLACEMENT 

FIG. 3.6 LOAD DISPLACEMENT CURVES FOR 
DIFFERENT CRACK LENGTHS. 
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In viev; of the above discussion, the J- ^ relationship 
can be studied, i,d.thout any loss of information, through the 
parameter (j/ 6^) representing slope of J vs lines in 
f; :'s, 3,4 and 3.5. 

For a particular laminate -with given crack 
configuration, the family of J curves can be reduced to 
a single curve of ( j/ 6^) plotted against crack length to 
specimen width ratio. 

For the purpose of analysis, j/ for the laminate 
was obtained with its primary stiffness matrix. After a ply 
failure, (l/ values were again obtained with the 
modified stiffness matrix. T.is process was continued till 
the laminate failed. 

Values of (j/ were obtained for quasi -iso tropic 
and crosspliod laminates having three different crack 
configurations, namely (i) SEN , (ii) DEN and (iii) centrally 
located crack, j/6^ values are given in Tables 3*2 to 3*7 
and shown graphically in Fi-~s. 3.7 to 3rl2. values 

have been shown for the laminate when all the plies are 
intact and also for the laminates after initial failure, i.e. 
when plies have failed in ascending order of their critidal 
displacement. These curves are used to predict the overall 
laminate behaviour. It may be noticed ( Tigs, 3*2 and 3,3) 
that for SEN specimens F vs, ^ curves for a/w ratios of 
0.5, 0.6 and 0,7 are quite close to each other for quasi- 
isotropic as well as crossplied laminates. This is in good 
agreement with experimental observations that J-integral curve 
is independent of crack length when a/w .>0,5 [18,20,22 1. 
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TABLE _3.2 
Single N’otched 

Quasi-iso tropic Laminate 

. .w . . .1 

, ALL PLIES I 90 PLIES , 90 and + 45 ° 

f ■ I failed plies faiied 

— - — .. ^ , 1 — I — - — 


O'.l 

; 9.78 

1 8.48 : 

6.10 ; 

0,2 

1 32.15 

' 27.82 

19.65 1 

0.3 

49 * 88 

43. 27 

29.55 i 

0.4 

■ 73. 25 

j 

63.62 

46.22 ! 

0,5 ! 

87.55 

: 76.22 

56.32 

0 , 6 j 

99.98 

f 

i 87. 20 

65.42 

0,7 ; 

: 102.0 

89.02 ’ 

67.48 

i 



TABLE 3.3 




Single Edge Notched 




Crossplied Laminate 


i 

a 1 


j/fi? in hj/m^-mm^ | 

I 

W i 

1 

- 1 

ALL : 

PLIES INTACT 

90 PLIES FAILED! 

0.1 


21.95 

18.- 2 

0. 2 


81.25 

1 71.5 

: 1 

1 0,3 


113.35 

94.35 

0.4 


154.60 

123.30 

0.5 


170.3 

131.6 

0.6 

1 

181 , 25 

! 136.15 

j 

0.7 

t 

[ 

178.55 

! 133.2 

I 



TABIE _3,L(. 

Double Jdge I'Totched 

luiisi-isotropic Laminate 


27 


2a 


tJ/ in 


vv 

ALL BLDG’S 

INTACT 

90° PLIES 
FAILED 

90° and + 45 ° 1 
PLIES FAT LED S 

0.1 

6.83 

5.92 

4.28 1 

j 

0. 2 

14.27 

12.28 

8.67 ! 

0.3 

20.88 

17.95 

12.77 1 

. O.lf 

27.18 

23.41 

1 16.68 

1 0.5 

I i 

33.09 

28.42 

i 

20 . 27 

1 0.6 

39.35 ! 

34 . 06 ' 

' 24. 28 j 

' 0.7 

\ 

i 

46.71 

40.52 

j 29. 00 ' 

1 i 

L I. 


TABLE 3,5 

Double Edge Notched 
Grossplied Laminate 


2a j 1 / 6^ in kj/m^ - mm^ 


1 j 

1 

ALL PLIES 

INTACT 

- - - 

90° PLIES I 

FAILED 1 

- - - - - - j 

1 

0.1 

15.35 

12.51 1 

0,2 

f 

36.27 

32.24 

i 0,3 

5 l.lif 

44.70 

1 

1 0.4 i 

66.74 

1 58,61 i 

1 1 

I 0,5 

79.53 

) 

! 68 . 91 1 

( I 

; 0, 6 

j 94.27 

! 80.96 ! 

i 1 

} 0.7 

i 

5 - - . 

108 . 25 

! 91.04 1 

1 i 

t \ 

i J 



TABLE _3.6 

Centrally Located Crack 
^'uasi -isotropic La'ulnate 


28 


2a ; 

¥ ; 

J / in 

kj/m^-mm^ 

i 

t 

f 

1 

t 

ALL PLIES I 

IkTACT 1 

90° PLIES 
FAILED 

90° and + 45® i 
PLIES FAILED ! 

0.1 1 

1 

5.90 i 

j 

5» 12 j 

1 

1 

3.66 1 

C\J 

• 

o 

12. 31 ■ 

10,66 

7.47 1 

0.3 1 

18. 51 i 

! 16.02 

: 11. 27 • 

0,4 

; 25.51 

1 

22,08 

15.56 ! 

0.5 

1 32,98 

28, 54 

20. 25 ; 

i 

0.6 1 

42.59 

! 36,82 1 

26. 34 i 

‘ ! 

0.7 

54.21 

! 

i 

; 46, 88 

• ! 

33.91 ; 

i 

i 
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T ABIE >3,7 


Centrally Located Crack 
Crossplied Laminate 


2a 

w 


j/ 6^ in kJ/m^-mm^ 


0,1 


ALL PLIES 
IlITACT 


Ilf-. 05 


0,2 


0.3 


0.4 


0.5 


0.6 


0.7 


33.06 
48,01 
65.29 
80,48 
99.22 
117, 88 




i 

j 

i 


90° PLIES 
EAIIEID 

12,08 
30.28 
42,94 
58.02 
69. 91 
83.97 
96.34 


■~n 



i 

i 

I 

I 

I 


j 

J 

1 





(J/d^), KJ/m 



a /w 


FIG. 3.7 J/6^ 


FOR QUASI - ISOTROPIC LAMINATE 



(J/6^); KJ/m^ 


£ 50 
E 


QUASI -ISOTROPIC 
CENTRALLY LOCATED CRACK 
' I - All Plies Intact 
II - 90° Plies Fdfled 
III - 90° and - 45° Plies Failed 



2a /w 

FIG. 3.9 J/6^ FOR QUASI -ISOTROPIC LAMINATE 


{J/6^) ; KJ/m 





(J/6^) ; KJ/m^ 


140p 

CR0S5PLIED (DEN) 

120- I - Plies Intact 

, II “ 90° Plies Failed 

Vi 

I 100 - 



2a/w 


FIG. 3.11 J/6^ FOR CROSSPLIED LAMINATE. 


(J/6^): KJ/m 



2a/w 


FIG. 3.12 J/6^ FOR CROSSPLIED LAMINATE 
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The procedure adopted for evaluating J-integral of 
the laminate through the stepwise procedure can he explained 
as follows- . Figure 3.13 shows an idealized load displacement 
curve of a specimen* The path O^TC * represents the model in 
which no stress relaxation takes place when a ply fails. 

The? path 04DEFG represents the model assuming complete load 
rolaxation upon a ply failure. For the analysis hero, a 
more piactical case where a partial load relaxation takes 
place ( represented hy OAI-IIJK ) has hoen assumed. The 
slopes Cj^, C 2 end c^ of the load-displacement curve segments 
between ply failures depend only on the initial crack 
length and are independent of the extent of load relaxation. 
If denotes a load rolaxation factor at a ply failure it 

can bo defined as follows ( with reference to Fig. 3*13). 

- 

''1 AD 

'^2 ^ ~ W ~ ^3.4-) 

To obtain J-integral at a given displacement * 6 area 
under the load displacement curve is first obtained . The 
throe areas shox'/n shaded in Fig. 3*13 can be easily 
obtained as 

= I =1 

^ = I <=2 "1 '> 

A^ = (l-r)j_) (cj_-C 2 ) (3*5) 

So the expression for potential energy can bo written as 
U = A^+ A^ + A^ = I + (l-n^)(c^-c2)6^ 

(3.6) 




6i 6 62 

DISPLACEMENT 


FIG. 3.13 LOAD -DISPLACEMENT CURVE FOR A LAMIN>j 
WITH PLY FAILURE. I 
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Nov 'J-inbC'-^r. i.l is rrivon by J 


a u 

^ a , so that 


J 


1 


J^C 

da 


QlC, 


1 1 2 P 

= 1 -I -a-^i) 


(--i _h_2 ) 6, ( 6.«,) 


;j a 


d a 

1 


(3.7) 


Defining kj_ = ^ | ^ and = - i 


da 


2 2 da 

2\ - . 


whoro kj^ and k2. arvi the values of (J/ ^ ) in these SG^Bionts, 
J= +2 (1 _Tl^)(k^ _.k2)&^(0 

= k2d^+ (k^ -k2)d^ +2(1 -il^)(k^ -k2)S^(& 

So if ’ ® • is in n segment then 
n 


(3.8) 


(?.9) 

Equation 3.9 has been used to obtain J- 6 plot for the 
laminates. 


3.3 COMPARISON WETII TIffi EXPERIMENTAL RESULT i 

I 

Predictions of for SEN specimens could be compared i 
with the experimental results which are available only for : 

j 

SEN specimens [20], Srinivasan [20] found that J vs d , 

curves for qua si -iso tropic as well as crossplied laminates t 
arc independent of crack length when ^ $>0,5. The 
experimental val. uc of for a laminate is obtained from this 
common J- 6 curve • In the present study, as has already been 
pointed out, J- 6 curves for SEN specimens, though change 
with ovack length, arc quite close to each other for 
a/w > 0,5. 
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In .}'f n‘’(\senfc study .1 values for the laminates 
has been obtained for three different models namely (i) 
when all t>n' plies fail simultaneously (ii) when plies fail 
successive Ly with no load relaxation in failed plies , and 
(iii) whe-n plies fail successively with complete load 
rclaxacion in failed plies* 

J vs curves for these throe models in the range 
0.5 ^ a/w <: 0,7 arc shown in Figs. 3*14 to 3,19 for (^pasi- 
isotropic and crossplied laminates. Values of obtained 
through different models are given in Tables 3,8 and 3,9 
alongwith the experimental values. It is observed that the 
experimental and theoretical curves are quite comparable. 
However, the model in vrhich all plies fail simultaneously 
predicts a higher value of J than the experimental value. 
This is expected because the off-axis plies do fail before 
the final fracture of the laminate. This model doosnot 
require a step-wise prediction method to bo applied and 
hence is simple to apply. 

The models in which plies fail successively predict 

lower values of than the experimental. But the model 

which assumes no load relaxation in failed plies predicts 

J- curve and J value closes to the experimental ones, 

c 

The difference between the predicted values of through 
the models with no load relaxation and with complete load 
relaxation is about 10^ and as may be expected former model 
predicts a higher value of Thus it may be concluded 

that it is basically the cxtensional stiffness matrix of the 




KJ/m2 


QUASI- ISOTROPIC (SEN) 
Successive Ply Failure 
No Load Relaxation 

0 Experimental 

1 a/w = 0.5 

II a/w=0.6 

III a/w=0.7 



p— w— *1 


I 1 


DISPLACEMENT, mm 

FIG. 3.15 J vs6 PLOT FOR QUASI - ISOTROPIC LAMINAT 






K J/ 



DISPLACEMENT, mm 

FIG. 3,16 J vs 6 PLOTS FOR QUASI - ISOTROPIC LAMINATE 



KJ/m2 



DISPLACEMENT , mm 

FIG. 3,18 J vs 6 PLOT FOR CROSSPLIED LAMINATE. 
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CROSSPLIED (SEN) 
Successive Ply Failure 

Complete Load Relaxation 

0 Experimental 

1 a/w = 0.5 

II a/w =0.6 

III a/w =0.7 



0.2 0.3 

DISPLACEMENT, mm 

FIG. 3.19 J vs 6 PLOTS FOR CROSSPLIED LAMINATE 
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0,5 

0,6 

0.7 


TABLE -3. 8 

Quasi -isotropic Liminate 


in 

Present 

kl/m^ 

work 


j "Expern” 

1 mental 

SirauTtaneous 
ply failure ; 

Successive ply failure 

no load Complete load 

relaxation ' relaxation 

! 

1 , 

i [20] 

! 

i 

If 9. 25 

" 1 
35.75 I 

31.95 

1 i 

1 

1 

i 

56# 24 

41. 30 ; 

37.10 

46.5 

57.37 

42# 46 ; 

1 

38. 25 1 



TABLE _3*9 
Crosspiied Laminate 



J'o i'’ 

kJ/m^ 




Present 

work 


Experim 

mental 

a/w 

Simultaneous i Successive 

1. . ... .. . 

ply failure 

1 ! 

0 

C\J 

1 — 1 


ply failure 

no load 

complete load 

! 

i 

|relaxn,tion 

1 

i 

relaxation 

1 

! 

r 

I 0,5 

I 

- - 1 

23.95 

20.58 

18.76 


i 0.6 

25.49 

21.56 

19.44 1 

22.65 j 

1 

O 

• 

25.11 

21. 16 

19.03 

' 

i 


4 
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lamin.-ice th, ),t controls its fracture toughness and stacking 
sequence mtiy not influence the fracture toughness signifi- 
cantly, The lower vaLuo of J predicted by stop-wise 
procedure may be attributed to the fact that in these models, 
Efjji and of the failed plies have been assumed to be 
ssoro and so in effect those layers carry no additional load 
in transverse diroction. However, in a practical situation 
these layers do carry additional transverse load away from 
tho dominant crack. 

It may be pointed out that the final fracture of the 
laminate was assumed to occur at the experimentally obtained 
critical displacement of the laminate and not at the critical 
displacement of the 0° plies. Although this limits tho 
applicability of the step-wise procedure , it is a necessary 
assumption in view of the experimental observations. 

Actually this is a consequence of a change in the critical 
displacement of 0° plius due to prior failure of off-axis 
plies and the resulting stress concentrations. 



CHAPTER IV 


-COWCLdSIONS AMD CUGC-ESTION FOR FUTURE WDRK 

J-integral f or quasi -iso tropic (0/+ 1^5 /qo)^^ and 
crossplied (0/90)j^g laminates has been determined by finite 
element technique, 

A stepwise analysis procedure was employed to obtain 
U^, After a ply failure stiffness matrix of the failed ply 
was modified for further analysis. Models with no load 
relaxation and complete load relaxation in failed plies 
were investigated. 

Based on the results presented in Chapter III the 
following conclusions are made, 

(i) U- ^ curve obtained for O.^ < a/w ^ 0,7 for SEN 
specimens are quite close to each other. This compares 
favourably with experimental finding of an unique 

1-6 curve for a/w 0,5. 

(ii) J values obtained by all the three models are in 

0 

good agreement with experimental result. The model in 
which all plies fail simultaneously predicts a higher 
value of J than experimental value. Models in which plies 

v> 

fail successively predict lower values of • The model 
in which no load relaxation in failed ply is assumed, ,fits 
the experimental results better than tho other two models. 



(ill) Difference between J values for the two models with 
plies failing successively is about lOf^, So that the effect 
of load relaxation in failed plies on J is marginal- 
This leads to the conclusion that extensional stiffness 
matrix of tho laminate controls the fracture toughness and 
stacking sequence may not influence it significantly, 

4. 1 SUGGESTIONS FOR FUTURE WORK 

(i) Experimental investigation for DEN and centrally 
located crack specimens should be conducted to find out 

I -integral for those specimens. 

(ii) A three dimensional analysis of laminates may 
be carried out to understand tho effect of interlaminar 
stresses and stacking sequence on J of the laminate. 

(iii) A failure criterion to predict ply failure in a 
cracked specimen should be established. 

Civ) J-intcgral should be obtained for laminates 
with off axis plies only, o.g. , angle ply laminates (+ e) 
and models presented in tho present study should be 
compared. 
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